Electronic correlations in graphite and carbon nanotubes from Auger spectroscopy 
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We have determined the screened on-site Coulomb repulsion in graphite and single wall carbon 
nanotubes by measuring their Auger spectra and performing a new theoretical analysis based on 
an extended Cini-Sawatzky approach where only one fit parameter is employed. The experimental 
lineshape is very well reproduced by the theory and this allows to determine the value of the screened 
on-site repulsion between 2p states, which is found to be 2.1 eV in graphite and 4.6 eV in nanotubes. 
The latter is robust by varying the nanotube radius from 1 to 2 nm. 
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Carbon nanostructures continue to be an intense field 
of both fundamental and applied research because of 
the recent discoveries of several of their unusual physi- 
cal properties. Among these one can recall (i) the ob- 
servation of the anomalous integer quantum Hall effect 
in planar grapheneP, 0] (H) the measurement of su- 
perconductivity at 11.5 K in Ca intercalated graphite 
and (in) intrinsic superconductivity in multi-wallQ and 
ultra-small Q carbon nanotubes at temperatures of 12 
and 15 K respectively. In the light of these unprece- 
dented properties and related new physics, the study 
and the quantitative estimate of electronic correlations in 
these carbon nanostructures are of paramount fundamen- 
tal importance. In fact, in one-dimensional conductors, 
like metallic nanotubes, the electronic interactions have 
a dramatic impact on their electronic properties, giving 
rise to the so-called Luttinger liquid behavior. This man- 
ifests in the power-law dependence of observables such as 
the tunneling density of states (DOS), of which suppres- 
sion at low energies has been observed in conductance 
measurements[l,[3- More importantly the accurate esti- 
mate of the screened Coulomb repulsion is a challenging 
problem that should be dealt within any theoretical study 
aiming at addressing the question of superconductivity. 

Auger electron spectroscopy is a powerful experimen- 
tal tool which permits the characterization of the effective 
interaction between electrons in solids. In particular the 
Auger lineshape is proportional to the 2-particle interact- 
ing DOS as a consequence of two valence holes creation on 
the same lattice site caused by the X-ray photoemission 
of a deep core electron. Several attempts have been made 
to interpret the Auger spectra of amorphous graphite [1| 
and highly oriented pyrolitic graphite (H0PG)[9[ but a 
satisfactory description is still to come. Moreover only 
few experimental data on single wall carbon nanotubes 
(SWCNTs) Auger lineshape are available Further- 



more, no theoretical effort introducing Coulomb repul- 
sion in SWCNTs has been attempted so far. 

In this paper we report on the comparative study of 
the Auger spectra of HOPG and SWCNTs. Through a 
new theoretical analysis of the Auger experimental data, 
we provide an accurate estimate of the on-site screened 
repulsion in both carbon structures. The access to this 
quantity is key in realistic local density approximation 
(LDA) -\-U calculations and in any low-energy interacting 
theory of the honeycomb lattice where only tt Dirac-like 
electrons are considered. 

SWNCTs were synthesized by ablating a CoNi-doped 
graphite target, using a pulsed Nd:YAG laser in the su- 
perposed double pulse configuration [l^. Raman spec- 
troscopy indicated that the tube are single wall, charac- 
terized by a low degree of defects and with diameters in 
the range 1.2-1.3 nm. [l3|. This is consistent with trans- 
mission electron microscopy (TEM) observations that, 
though showing tubes aggregated in bundles of various 
dimension and twisting, allowed us to measure a tube di- 
ameter of 1.2± O.lnm through a statistical analysis 11|. 
Moreover electron energy loss spectroscopy performed, 
by using the TEM apparatus, directly on SWCNTs bun- 
dles at the Co and Ni L ?, 3 edges did not detecte any 
traces of these catalysts. [l2| A droplet of the synthesis 
product was diluted in isopropyl alcohol and dispersed 
on a metallic surface. A freshly cleaved HOPG sample 
was used for measuring the core- valence-valence (KVV) 
Auger features. The Auger spectra were acquired using 
an Al Ka (1486.6 eV) monochromatic x-ray source with 
a resolution of about 1 eV. The obtained experimental 
spectra are shown in Fig. 1 after subtraction of secondary 
electron background. 

The Auger lineshape of solids can be calculated by us- 
ing the so-called 2-step approach, in which the photoe- 
mission and the Auger decay are considered as indepen- 
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dent events. In absence of significant electronic correla- 
tions, the computation of KVV Auger spectrum reduces 
to the self-convolution of the 1-particle valence DOS. If 
moderate or strong (compared to the bandwidth) on-site 
repulsion is present, the lineshape can be calculated by 
means of the Ci nill3l |-Sawatzkv[14| approach. 
Following Ref.[13j], the Auger current J reads 



J ^ ^ ,a2 ,03 ,a4 ,cr-^ai ,a2 ,a3 ,04 ,cr (^) (1) 

Ql ,a2 ,a3 ,Q4 ,CT 



where Ui denote all the single-particle valence orbitals 
available in the solid, A is the so-called Auger matrix 
element given by 
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X (w|4<T,dfe^J^Ml3|44^b) (2) 



with k, c and Uk, cr^ denoting the Auger electron and core 
orbitals and spin respectively. D is the 2-particle inter- 
acting DOS 

' (3) 

where H is the interacting hamiltonian of the solid. Here 
we denote by \v) the hole- vacuum and by d-^^ the anni- 
hilation (creation) operator of a hole in spin-orbital i. D 
is obtained as usual from the anti-hermitian part of the 
2-particle Green's function Gai^a2,a3,ai,a{'^) which obeys 
the matrix Dyson jT6| equation 



= GW[l + c/,G(^')]- 



(4) 



where G^"-* is the noninteracting 2-hole Green's function 
and U is the matrix of screened on-site repulsion for va- 
lence states. The screened interaction differs from the 
bare atomic one, defined as 



U. 



Oil ,0:2 ,0:3 ,CK4,CT 



{v\dc,,^da^„\ — \d\^^d'^^^„\v). (5) 



The evaluation of U starting from the atomic value U'^ is 
generally a delicate task. In the following we discuss the 
phenomenological approach we have adopted to deter- 
mine this quantity. The Cini-Sawatzky approach works 
quite well in closed (or almost closed)-band systems like 
zinc and copper, where the ladder approximation lead- 
ing to Eq. (HI) provides an exact result. However, if the 
Fermi level crosses the middle of the conducting band, the 
computation of the Auger current becomes a remarkably 
challenging many-body problem, which usually cannot 
be solved by evaluating Green's functions [iHj . 

In the fight of this, the theoretical study of Auger spec- 
tra of HOPG and SWCNTs is indeed far from straight- 
forward because the a and tt bands are half filled. How- 
ever in these systems some special features (which are 
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FIG. 1: (Color online) Experimental KVV Auger spectra of 
HOPG graphite (bold red curve) and SWCNTs with average 
diameter of 1.3nm (black curve) 



discussed hereafter) allow the use of closed-band theory, 
with slight but crucial modifications. 

First we observe that the DOS is largely suppressed 
in the proximity of the Fermi level, so that screening is 
not very efficient. This implies a static renormalization 
of the bare interaction [/'' which must be used in the 
theory. Second, we recall that the bonding portion of 
the (Ts,p bands is separated by several eV from the anti- 
bonding part located above the Fermi level. As long as 
such a separation is larger than the effective interaction, 
one can treat the band as if it was closed, thus justifying 
the approach reported by Cini[l3|. where no structural 
modification is needed for the interacting Green's func- 
tion in Eq. (|4]). However, the situation is different for 
the TT band, where the bonding and antibonding portions 
are separated by a very small region with a small DOS. 
Here Cini's approach can not be used without appro- 
priate modifications. In this case the contribution to the 
Auger spectrum originating from tt and mixed it — a holes 
would be strongly influenced by open-band effects. It is 
also expected that such a region should reveal the princi- 
pal differences between the spectra of HOPG and SWC- 
NTs. In fact screening and excitonic effects 17[ and Lut- 
tinger liquid properties in SWCNTs are expected to lead 
to a quite different behavior of electrons in proximity of 
the Fermi level due to the different dimensionality. This 
conjecture seems to be confirmed by the experimental 
data. Indeed, the tt and mixed tt — ct portion of spectrum 
(i.e. ljj > 250 eV) show clear differences between HOPG 
and SWCNTs, while in the region (i.e. lo < 250 eV) 
the two spectra are quite similar. In particular for 250 eV 
< w < 280 eV the lineshape of SWCNTs is narrower with 
vanishing and much weaker intensity in proximity of the 
Fermi level, as compared to the one for graphite. This 
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FIG. 2: 1-particle partial DOS of (10,10) SWCNT (diameter 
close to 1.3 nm) obtained by the tight binding method of 
Refl2^. The inset shows the same quantity for graphite, taken 
from Ref@. The Fermi level corresponds to zero-energy and 
the antibonding part is not displayed. 



fits well with a scenario where the screening properties 
of TT electrons are less efficient in SWCNTs. 

Within the closed-band theory, the Auger spectrum 
is obtained by taking the Auger matrix elements and 
the on-site interactions from atomic calculations, which 
neglect solid state effects. On this basis, one intro- 
duces the static screening operated by the closed-band 
system simply by rescaling all the F^^\i,j) Slater in- 
tegrals that enter the bare U'', such that F'^^^{i,j) 
F'-°\i,j) -W. W can be taken as the unique free fit- 
ting parameter of the theory. Alternatively W can be 
also estimated within the Random Phase Approximation 
or ab initio methods [is']. The only ingredient which ac- 
counts that the Auger holes are in the solid is the non- 
interacting 1-particle DOS p*^°^(a;). Its self-convolution 
D^^''{lu) — J de p'^'^\e) p'''^\lo — e) and the corresponding 
Hilbert transform build the noninteracting G'^ entering 
Eq. dH). 

Cini's approach should in principle be completed by 
introducing the effect of off-site interaction. Experiments 
on Au [l^ showed that there is a shift of 1.2 eV between 
the profile predicted by the above theory and experiment. 
The shift is 2.4 eV in the case of Ag ^.20J. This was 
explained in terms of the off-site interaction. In the two- 
hole resonance there is an important amplitude that the 
holes sit on neighboring sites, and including the nearest- 
neighbor interaction into the theory yields an almost rigid 
shift close to the experimental one[2lj. 

In the following we will phenomenologically consider 
the open-band effects by introducing orbital-dependent 
form factors fai,a2,a3,a4- This must be introduced to 
correct all the quantities measuring local properties ex- 



pressed by {v\daida2\0\dl^^dl^Jv) where O is a local 
observable. Therefore the effective on-site repulsions 
Uai,a2,a3,ai,a (whcrc has been already rescaled by 
W) and the matrix elements Aa-^^a2-aa,a4,.tT will be cor- 
rected by a comnioii multiplying factor fai,a2,a3,a4- 

In 

our case the states are a^, da;, tiy, tt. The form factor 
/ takes into account that the 2s states of carbon behave 
as if they were atomic, while the 2p ones are delocalized 
in the lattice. The latter can use only 1/2 of the total 
CTp and TT states to form occupied localized states because 
the p-bands are half-filled. Therefore we have three in- 
dependent / factors corresponding to having (i) four as 
orbitals , (ii) two as and two ax.y^T^ orbitals, (Hi) four 
ax.y,Tr orbitals in the quartet {ai, q;2, as, Q!4}. Accord- 
ing to the above discussion the three independent form 
factors are estimated to be fssss ~ 1, fsspp ~ 1/2 and 
fpppp ~ 1/4. We will show that this choice works quite 
well in the case of HOPG, while we need fpppp ~ 1/2 
to reproduce the Auger spectrum of nanotubes. Indeed 
in nanotubes the geometry constrains the holes and this 
could be the reason for a larger fpppp than in grap hite. 
It is worthwhile to note that the analysis of Ref.[21] does 
not apply to p holes and in fact no shift is seen in this 
case (the pairs presumably extend further than a near- 
est neighbor distance). A shift could be present in the 
KLiLi case, but we cannot tell since there is a single 
peak there. 

We proceed by evaluating the noninteracting 1-particle 
DOS p*^"^ for each kind of valence state. In the case 
of HOPG, we use the DOS from Ref^S which is taken 
from experiments. For SWCNTs we performed a tight 
binding calculation [23| including both 2s and 2p or- 
bitals, but neglecting overlap integrals for simplicity. 
The result for a typical (10,10) armchair nanotube 
with diameter close to 1.3 nm is shown in Fig[2] to- 
gether with the DOS of HOPG. For the Auger ma- 
trix elements, we used the (spin-independent) values 



As 



0.8, ^ 



0.5, A 



1.0 which are ob- 



tained by atomic calculations [g| and hence apply to both 
graphite and carbon nanotubes. The bare (atomic) 
on-site Coulomb repulsions are obtained by appropri- 
ate combinations of the Slater integrals F^^'^'>{i,j) and 
G^'^Hhj)^ found in the literature [lil. The indepen- 



dent bare interactions are (in eV) U^sssi — 15-5, C/^, 



1.5. 



15.0, U^xxsi 

Tjb — Tjb 

SXSX^ ' XTTX-K^ 



Tjb (- 
TTTrTTTr J. V- 

XTTTTxi 



U'' ) 



14.6, U', 



sxsxl 
I 

13.9, U^^^^^ = 0.8, C/^-^^ =^1.9, 
— 13.1. As discussed above, these 
values must be corrected by solid state effects. This is 
done by subtracting the screening constant W from the 
P^'^Hhj) Slater integrals and multiplying the resulting 
U and A matrix elements by the / factors {W being the 
the only fitting parameter of our approach). 

The theoretical spectra of HOPG and SWNCTs were 
computed by solving a 16 x 16 matrix problem for a =i 
and a 6 X 6 problem for a =|, as shown in Eq. (|H)[l^. 
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FIG. 3: (Color online) (a) Theoretical lineshape [computed 
from Eq. HJ] of KVV Auger spectrum for HOPG (bold red) 
and for SWCNTs (black) curve; (b) Diagonal contributions of 
the interacting DOS for HOPG (bold red) and for SWCNTs 
(black) where the two valence holes have the same symme- 
try. The DcTpCTpcrpcrp Contribution is understood as the sum 
^CT^<T^<T^<T^ + DtTyCTyCTyay Thc a;-axis displays kinetic energy, 
obtained by shifting the position of the Fermi level in Fig. 2 
of 284.6 eV, which is the binding energy of Is core hole. 



The final result is plotted in Figl3^, where the best fit- 
tings yielded the respective values Whopg = 6.0 eV and 
W^swcNT = 5.5 eV for HOPG and SWCNTs. The agree- 
ment between theory and experiment is quite good, and 
is particularly satisfactory for graphite. 

These values permit the determination of the most rel- 
evant parameter of our model, which is the screened on- 
site repulsion between the 2p states. Thus, the best 
fitting for W yields Upppp = 2.1 eV for HOPG and 
Upppp = 4.6 eV for SWCNTs. This result gives rise to 
the lack of features close to the Fermi level for SWCNTs, 
making the Auger spectrum more symmetric and nar- 
rower than that of HOPG. This is understood by looking 
at FiglHb, which shows the diagonal contributions of the 
interacting DOS according to Eq. ([3]), where the valence 
holes were taken in the same state. The off-diagonal con- 
tributions are not shown for the sake of clarity but are 
essential to reproduce the experimental spectra. 

Concerning the lineshapes, the most striking feature 
is the narrow structure at 240 eV, which also appears 



as a shoulder in the spectrum reported by Houston et 
al.[§|. This peak was assigned to a plasmon rephca of 
the main structure at 265 eV produced by a plasmon 
with an energy lOp — Tl eV. Conversely we interpret 
the narrow structure as a quasi-two-hole resonance pro- 
duced by two Og Auger holes. This is consistent with 
the predicted values of the screened on-site repulsion be- 
tween Gs holes, which are Ussss = 9.5 eV and 10.0 eV 
for HOPG and SWCNTs respectively The noninteract- 
ing D^sL has a maximum at Sss = 252 eV (graphite) 
and 251 eV (nanotube) and therefore a narrow structure 
around Sss — Ussss ~ 241 eV in the interacting Dgsss is 
correctly expected. Since the CTs-bandwidth is of ^ 20 
eV a full splitofF two-hole resonance cannot happen, but 
a strongly distorted band-like behavior occurs (see Fig. 
lb of RefHsl). It is worth noting that Auger spectrum 
from a sample consisting of SWCNTs with average di- 
ameter of 2nm does not show significant changes with 
respect to that reported in FigUJ Moreover, by perform- 
ing a similar theoretical analysis on a (20,20) SWCNT 
no substantial changes can be found for the values of the 
correlation interaction. This means that the values we 
obtain for the correlation in SWCNTs have a very small 
dependence on the nanotubes diameter. 

In conclusion the lineshape of the Auger spectra for 
HOPG and SWCNTs have been interpreted in terms of 
a new theoretical approach using a single fitting param- 
eter. The Upppp Coulomb repulsion results doubled pass- 
ing from HOPG to SWCNTs. This explains the sizeable 
shift of the Auger feature at high kinetic energy measured 
for SWNTs, as compared to HOPG. Finally we point out 
that the increase of the Upppp value is consistent with the 
theoretical prediction [2^ of the enhancement of the su- 
perconductive critical temperature observed recently in 
carbon nanotubes. 
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